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ABSTRACT 

The heating associated with the deposition of 7-rays in an accretion disk 
is proposed as a mechanism to facilitate the transformation of a low mass X- 
ray binary to the radio millisecond pulsar phase. The 7-ray emission produced 
in the outer gap accelerator in the pulsar magnetosphere likely irradiates the 
surrounding disk, resulting in its heating and to the possible escape of matter from 
the system. We apply the model to PSR J1023+0038, which has recently been 
discovered as a newly born rotation powered millisecond pulsar. The predicted 
7-ray luminosity ~ 6 x 10 34 erg s -1 can be sufficient to explain the disappearance 
of the truncated disk existing during the 8 month~ 2 yr period prior to the 2002 
observations of J1023+0038 and the energy input required for the anomalously 
bright optical emission of its companion star. 

Subject headings: binaries: close — magnetic fields — pulsars: general — pulsars: 
individual (PSR J1023+0038) — stars: neutron 
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1. Introduction 

The discovery of pulsed radio emission at 1.69 ms from the source J102347.67+003841.2 
(see Archibald et al. 2009), offers a unique opportunity to study the evolution of a system 
from the low mass X-ray binary (LMXB) phase to the recycled millisecond pulsar (MSP) 
phase. The system has an orbital period of 4.75 hrs (see Woudt, Warner, & Pretorius 2004) 
and is distinguished from other systems by a change in optical appearance from a spectrum 
characterized by broad double peaked hydrogen and helium emission lines (see Bond et al. 
2002) to a late type absorption spectrum (see Thorstensen & Armstrong 2005). Further 
analysis based on SDSS spectra by Wang et al. (2009) provide evidence for the presence of 
a truncated accretion disk in the early observations in 2001 and its possible absence in 2002. 
The recent detection of X-ray pulsations by Archibald et al. (2010) suggests that pulsar 
high-energy emission activity was already present in 2008. 

The description of the transition from the accretion powered to the rotation powered 
phase is not understood due to the complexities in the description of the interaction between 
the magnetosphere of a neutron star (NS) and its accretion disk. Simple arguments have 
been discussed in early work by Campana et al. (1998). Although theoretical developments 
have not progressed significantly, recent observations of the optical counterparts of NS soft 
X-ray transients in the quiescent state have provided indirect evidence for additional heat- 
ing of the companion star. As examples, orbital modulations in the optical emission from 
J102347.68+003841.2 (Thorstensen & Armstrong 2005), SAX J1808.4-3658 (Burderi et al. 
2003; Deloye et al. 2008), XTE 1814-338 (D'Avanzo et al. 2009) and IGR J00291+5934 
(Jonker, Torres & Steeghs 2008) have been found in their respective quiescent states. Specif- 
ically, the amplitude of the modulation cannot be explained by irradiation associated with 
the X-ray emission from the disk or NS due to insufficient luminosity, indicating the need 
for the operation of an additional heating source. To provide an explanation of the orbital 
modulation of the optical emission, pulsar wind models have been suggested in which heat- 
ing of the donor is due to the effect of a relativistic pulsar wind (Burderi et al. 2003). In 
addition, the interaction between the pulsar and its companion has been suggested as leading 
to the formation of an intrabinary shock (see Campana et al. 1998) for the generation of 
non-thermal X-ray radiation observed in J0024-7204W in the globular cluster 47 Tuc (see 
Bogdanov, Grindlay & Berg 2005), in PSR J1740-5340 in NGC 6397 (Grindlay, et al. 2002) 
and in PSR J1023+0038 (Archibald et al. 2010). Finally, the disappearance of the truncated 
disk in J102347.67+003841.2 may have been a direct result of the activation of the pulsar. 

In order to facilitate this transition, the accretion of matter onto the NS must decrease 
rapidly, for otherwise, the NS will spin down (e.g., see Jeffrey 1986). In contrast to the case 
of radio MSPs in long orbital period systems where a red giant companion can detach from 
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its Roche lobe as its envelope is exhausted, a sudden decrease in the mass transfer from the 
donor is not generally expected for the donors in short period binary systems. For these 
systems, the disk must be truncated sufficiently such that the mass flow is prevented from 
entering into the NS magnetosphere. A mechanism possibly leading to a sudden decrease in 
the accretion flow is the operation of the so called "propeller" effect (Campana et al. 1998; 
Romanova et al. 2009), in which the Alfven radius exceeds the corotation radius, during the 
quiescent state of a X-ray transient LMXB phase. In this regime, however, the prevention 
of matter accreting onto the NS is uncertain. That is, matter can either be swept away 
or accrete intermittently onto the NS. Ablation due to the action of a pulsar wind may 
also facilitate the removal of the disk around the pulsar (Wang et al. 2009), however, an 
important theoretical uncertainty is the description by which a particle kinetic dominated 
flow can be formed near the pulsar (Kirk & Skjaeraasen 2003; Arons 2008). 

In this Letter, we explore the possibility that irradiation by 7-ray emission produced 
in the outer gap accelerator within a pulsar magnetosphere can facilitate dissolution of the 
disk. In §2, we briefly describe the theoretical concepts of the model and apply the theory 
to J102347.67+003841.2 in §3. In the last section, we discuss the implications of our model 
and conclude. 



2. Theoretical Model 

The activation of the radio pulsar phase of a NS in a short period LMXB is likely to 
occur during the quiescent state (see Campana et al. 1998) when the mass accretion rates 
are lower than the time averaged value, typically M ~ ]_q14-16 g^/ g (Heinke et al. 2009). 
In such a phase, the magnetosphere can extend beyond the light cylinder radius. If the 
pulsar magnetosphere is sufficiently clear of matter, 7-ray photons with an energy ~ GeV 
can be produced in an outer gap accelerator (Cheng, Ho & Ruderman 1986a, 1986b) via 
the curvature radiation process of high-energy particles. The emitted 7-rays can irradiate 
the accretion disk to promote the escape of disk material from the system (c.f. Figured]). 
Sufficient heating of the disk may directly lead to the escape of matter or to an increase in 
its vertical scale height. In this case, the disk would subtend a larger solid angle as seen from 
the NS, thereby facilitating the acceleration of disk material associated with the radiation 
pressure of the magnetic dipole radiation. In the following, we describe the process. 
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2.1. 7 _ra y emissions from the outer gap 

In the outer gap model for the rotation powered pulsars, particles can be accelerated by 
the electric field along the magnetic field lines in the region where the local charge density 
deviates from Goldreich- Julian charge density. The typical strength of the accelerating field 
in the gap is expressed as (Zhang & Cheng 1997) 

where V a = B S R 3 S /Rf c is the electrical potential drop on the polar cap, P is the rotation 
period, B s is the stellar magnetic field, s is the curvature radius of the magnetic field line, 
Ric = cP/2ir is the light cylinder radius, and R s is the stellar radius. In addition, / is 
the fractional gap thickness, which is defined by the ratio of the gap thickness to R\ c . By 
assuming force balance between the acceleration and radiation reaction of the curvature 
radiation process, the typical Lorentz factor of the particles in the gap is 

r = (f^ii) V4 ~ 1.6 x i^^p-TbT^rT, (2) 

where P_ 3 = (P/l ms), B 8 = (B s /10 s G), s 1 = (s/R lc ), and R 6 = (P s /10 6 cm). The 
corresponding photon energy of the curvature radiation is 

E, = ^ ~ 25 / 3 / 2 P~3 7/ V V 1/4 ^6 /4 GeV. (3) 

The 7-rays emitted in the gap will be converted into electron/positron pairs via the photon- 
photon pair-creation process with X-ray photons from either the NS surface or the accretion 
disk. Because the created pairs tend to screen the electric field in the gap, the gap thickness 
is controlled by the photon-photon pair-creation process. Using the pair-creation condition 
that E^E X ~ (m e c 2 ) 2 , / is estimated as 

/ ~ 0.22PZ{ 6 5- 1/2 *l /e X 3/2 Erf 3 , (4) 

where £0.1 is the typical X-ray photon energy in units of 100 eV. The luminosity of the 7-ray 
emission from the outer gap is estimated as 

L 7 ~ fL sd ~ 4.2 x 10 33 P„- 3 1/2 5 8 1/2 S ; /2 ^ /2 P :i 2 erg/s, (5) 

where L sd = 4(27r) 4 B 2 R 6 / '6c 3 P 4 is the pulsar spin down power. 

The rotation period of the newly switched on MSP is related to the history of the 
accretion onto the NS. In particular, its rotation period in the accretion stage may be related 
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to the equilibrium spin period, which is obtained by equating the co-rotation radius, r co = 
(GMP 2 /47r 2 ) 1/3 to the magnetospheric radius, r M ~ 2.4 x 10 6 B* /7 Rl V7 M~ 2/7 M~l /7 cm. 
Here, M is the NS mass, M 15 is the time averaged accretion rate in units of 10 15 g/s and 
Mx.4 is the NS mass in units 1.4 M . However, the spin period may not be equal to its 
equilibrium spin period since the pressure build up associated with the accumulation of disk 
matter at the magnetospheric boundary can push the magnetospheric boundary inside the 
magnetospheric radius r« (Romanova et al. 2009), resulting in a higher spin period than its 
equilibrium value. Because of this uncertainty, we introduce a parameter, (3, as r co = (3tm, 
which yields the rotation period, 

P e ,_ 3 = 1.7kB 6 8 /7 . (6) 

where k = (3 3 / 2 Rf /7 M~ 3/v M~ 4 5/r . The spin period is illustrated in Fig. [2] clS Sl function of 
the stellar magnetic field with k — 1 (solid line), 0.4 (dashed line) and 0.2 (dashed-dotted 
line). For example, with (3 = 0.5 and Rq = M1.4 = 1, k — 1, 0.4 and 0.2 corresponding to 
M 15 ~ 0.1, ~ 0.7 and ~ 4. In addition, we also plot the observed spin periods and inferred 
magnetic fields of the radio MSP population in the Galactic field. Only MSPs in binaries 
with orbital periods less than about 10 days are included since those systems with longer 
periods may have had red giant companions. In this case, the dissolution of the disk is due 
to the cessation of mass transfer as the envelope of the giant star is exhausted leading to the 
shrinkage of the red giant within its Roche lobe. 



2.2. 7-ray irradiation of the disk 

The magnetospheric 7-rays irradiating the surrounding disk may be absorbed via the 
so called pair-creation process in a Coulomb field by nuclei in the accretion disk. Further 
absorption will occur as the relativistic pairs created with a Lorentz factor T ~ 10 3 will 
transfer their energy and momentum to the disk matter via the Coulomb scattering process. 

The cross section of the above pair-creation process is given by Lang (1999) as 



<x(£ 7 ) = 3.5 x 10" 3 ZV r 
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where ot is the Thomson cross section, Z is the atomic number, a is the fine structure 
constant, and m e c 2 is electron rest mass energy. All 7-rays irradiating the disk can be 
absorbed if the column density of the accretion disk exceeds a critical value given by 



m p h 



a I 



h 2 
60y g/cm , 
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where I is the propagation length of the 7-rays in the disk, h is the thickness of the disk, 
Z 2 = 3 appropriate for a solar abundance and m p is the mass of a proton. It can be seen 
that the above column density is smaller than the typical column density of the accretion 
disk. For the standard gas pressure supported disk model (Frank, King & Raine, 2002), the 
column density is 

E(Ric) ~ 2.1 x lO^MTMl^PT g/cm 2 , (9) 

where «o.i is the viscosity parameter in units of 0.1. We find that the required column 
density OH]) is much smaller than that of the standard disk model for viscosity parameters 
less than unity. 

An estimate of a critical accretion rate below which the radiation associated with 7- 
rays from the outer gap can energetically eject the matter can be obtained by equating 
L 7 = GMM/2R lc with equation ©, given by 

M c = 2.1 x 10 u P^Bl /2 sl /2 Rl /2 Eo,lMi l g/s. (10) 

Applying the critical accretion rate ( ITUj) and (3 = 0.5 to equations §6§ and using equa- 
tion the predicted 7-ray luminosity and the surface magnetic field are illustrated as a 
function of equilibrium period and the X-ray photon energy in Fig. [3j The different lines 
represent the results for the different X-ray photon energies. For a pulsar with P_3 = 5 and 
E i = 0.5, for example, the predicted 7-ray luminosity and the surface magnetic field are 
L 7 ~ 5 x 10 34 erg s _1 and Bg — 41, respectively. Note that if the true accretion rate is less 
than the critical accretion rate, the predicted magnetic field and 7-ray luminosity for the 
specific set of the rotation period and the X-ray photon energy provide smaller value than 
those associated with the critical rate. 



3. Application to PSR J1023+0038 

As an application of our model, we estimate the 7-ray luminosity which inhibits accretion 
into the magnetosphere, once the rotation powered pulsar is activated. In Fig. [2j the observed 
characteristics (filled box) for PSR J1023+0038 are indicated. For the X-ray emission, Homer 
et al. (2006) found that the X-ray emission in 2004 observations are well fit by a NS 
atmosphere model (Ex ~ 30 eV) plus power law component. However, we choose the X-ray 
emission from the NS stellar surface rather than the X-ray emission from the pair-creation 
originating from either the pulsar magnetosphere or a shock where material overflowing the 
companion interacts with the pulsar wind. If the X-ray emission is magnetospheric in origin, 
the synchrotron emission from the secondary pairs produced outside the main acceleration 
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region is possible. In such a case, because the high-energy photons are propagating to convex 
side of the magnetic field lines, the non-thermal X-ray photons are unlikely to illuminate the 
gap region and, therefore, collisions with the 7-rays inside the gap are not expected. On the 
other hand, if the non-thermal X-rays originate from the shock region, the photon-number 
density in the pulsar magnetosphere is much smaller than that associated with the thermal 
X-rays from the NS surface, indicating the photon-photon pair-creation process in the gap 
mainly occurs with the thermal X-rays. 

As Fig. [2] reveals, the parameter k < 0.4 in the relation ([H]) can be consistent with 
the observed rotation period and upper limit of the inferred surface magnetic field. This 
may indicate that the averaged accretion rate is order of M ~ 10 15 g/s unless /3 is much 
smaller than unity. On the other hand, the predicted 7-ray luminosity (JSJ) and the critical 
accretion rate ffTUl) give L 1 ~ 6 x 10 34 erg s _1 and M ~ 5 x 10 15 g/s, respectively, with 
B$ = 3 and Eq.i = 0.3 and s\ = R§ = M1.4 = 1, indicating the 7-ray irradiation on the disk 
would be sufficient to inhibit the accretion flow into the magnetosphere. We note that the 
predicted luminosity L 1 ~ 6 x 10 34 erg s _1 is consistent with the required luminosity from an 
irradiating source to explain the heating of the companion star (Thorstensen & Armstrong 



As argued by Wang et al. (2009), there is evidence that a truncated disk existed in 
SDSS J102347.6+00381 during the 8 month~ 2 year period prior to the 2002 observations. 
Based on a standard disk model a disk extending from 10 9 cm to 5.7 x 10 10 cm with a total 
disk mass of ~ 1.7 x 10 23 g was inferred. Here, we examine whether the purported 7-ray 
irradiation can remove the disk within a time scale of less than 2 yr. 

Wang et al. (2009) estimate a disk column density of ~ 36(a/0.1) _4 / 5 M 1 7 / 10 r ] ~ 3//4 , 
which is higher than the required density ([H]) for 7-ray absorption because the propagation 
distance in the disk is much longer than the disk thickness. Here, we examine the fraction 
of the 7-ray power which can irradiate the disk at r = 10 9 ~ 10 10 cm. If we apply (1) 
a pure dipole magnetic field, (2) the geometry of the outer gap accelerator, which extends 
from the so called Goldreich- Julian null charge surface, where the product Q ■ B — 0, to 
the light cylinder, and (3) a standard disk model, implying a disk thickness of h ~ 10 s cm 
at r ~ 10 10 cm from the pulsar, we find that about 77 ~ 0.1 % of the outer gap region can 
illuminate the disk extending beyond r > 10 9 cm. Assuming the 7-ray irradiation can remove 
the matter at a rate of M = rjL^r/GM, the time scale of the life of the disk, r ~ M^/M, is 



2005). 




(11) 



which can be consistent with the life time inferred from the observation. 
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4. Discussion & Conclusion 

Since our scenario for describing the phase from the accreting MSP phase to the radio 
MSP phase invokes the 7-ray radiation from the outer gap, observations using the Fermi 
satellite provide important constraints on the 7-ray emission model. Since the spin down 
power of PSR J1023+0038 (L sd < 3 x 10 35 erg/s) is relatively high compared to 7-ray 
emitting pulsars detected by Fermi (Abdo et al. 2009) and its distance (d ~ 1 kpc) is 
similar to these Fermi sources, PSR J1023+0038 is a prime candidate for 7-ray emission. 

Recently Tarn et al. (2010) report the discovery of the 7-ray emissions in the direction 
of PSR J1023+0038 with a flux 5 x 10~ 12 erg curV 1 . Although the pulsation has not 
been confirmed in 7-rays by Tarn et al. (2010), its position is in coincidence with the radio 
position, which lies within the 68% confidence-level error circle of the 7-ray position. The 
7-ray flux can be estimated as _F 7 ~ L 1 /5VLd 2 ~ 10 -9 erg/cm 2 s, where L 7 = 6 x 10 34 erg s _1 , 
the typical solid angle SQ = 2, and the distance d — 1.3 kpc are assumed. The theoretical 
prediction is found to exceed the measured value, however, the 7-ray flux is dependent on 
the viewing geometry. For example, in the the outer gap model the 7-ray emission decreases 
with decreasing viewing angle as measured from the rotation axis. 

Given the constraints on the inclination of the orbital plane of PSR J1023+0038, ~ 34° 
to ~ 53°, it is likely that the line of sight will pass through edge of the 7-ray beam from the 
outer gap accelerating region, and that the observed 7-ray flux will be much smaller than 
_F 7 ~ 10~ 9 erg cm -2 s -1 . Because the 7-rays associated with edge of the beam will be emitted 
with an accelerating electric field smaller than the typical value described by equation (JTJ), 
the energy of the 7-rays will be smaller than that given by equation ([3]). Hence, the recent 
measurement of 7-ray radiation in the direction of PSR J1023+0038 can be reconciled with 
the proposed emission model provided that the line of sight passes through the edge of the 
7-ray beam. With the typical outer gap geometry, the 7-ray pulse profile is expected to be 
either a single peak or double peak with a small phase-separation for this viewing orientation 
(c.f. Hirotani 2008). 

The secondary pairs produced outside the gap region by the pair-creation process will 
emit non-thermal X-rays via the synchrotron radiation process with a typical energy of 
E n ,x ~ 2(r/10 3 ) 2 (i? sin 8/10 5 ) keV, where T is the typical Lorentz factor of the secondary 
pairs, B is the magnetic field near the light cylinder and 9 is the pitch angle. The optical 
depth of the pair-creation process of a 7-ray can be estimated as 7"x 7 ~ ^s^x-yRic ~ 10~ 2 
with n, % x ~ L S)X / '(AttRi c cE Sj x) corresponding to the X-ray photon density, where we adopted 
the observed X-ray luminosity L S) x ~ 10 31 erg s -1 and E s> x ~ 30 eV (Homer et al. 2006), In 
addition, <7x 7 ~ or/3 is the cross section of the photon-photon pair creation. The luminosity 
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of these non-thermal X-rays is estimated as 

L n , x ~ r X7 L 7 ~ 6 x 10 32 erg/s. (12) 

which is consistent with the observations L n ^x ~ 5 x 10 31 erg/s (c.f. Homer et al. 2006) 
assuming that we may be observing the edge of the high-energy emission region. The spectral 
shape below E n>x ~ 2 keV will be expressed by the photon index T n X ~ 1-5, which is near 
the observed value T n> x ~ 1.3 reported by Homer et al. (2006). 

To conclude, it has been shown that the irradiation by 7-rays from the outer gap can 
facilitate dissolution of an accretion disk within the assumption that the rotation powered 
pulsar is active during a sufficiently low luminosity level in the quiescent state of a soft X-ray 
transient system and the heating leads to the escape of disk material from the system. The 
model has been applied to the recently born MSP PSR J1023+0038, and we have estimated 
that the 7-ray luminosity is ~ 6 x 10 34 erg s _1 , which is consistent with the irradiated 
luminosity required to explain the heating of the companion star. 

Future 7-ray observations of the pulsed component are needed to confirm that J1023+3008 
is a 7-ray emitting pulsar and to place constraints on our theoretical model of the 7-ray emis- 
sion, its site of production, and its effectiveness in facilitating the removal of the accretion 
disk in this system. 

This work was supported in part by the Theoretical Institute for Advanced Research 
in Astrophysics (TIARA) operated under the Academia Sinica Institute of Astronomy & 
Astrophysics in Taipei, Taiwan and by NSF AST-0703950 to Northwestern University. 
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Fig. 1. — The schematic view of irradiation by the 7-rays from the outer gap (shadowed 
region) within a pulsar magnetosphere. The outer gap accelerator is extending between the 
null charge surface, where Q ■ B = 0, to the light cylinder. 
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Fig. 2. — Plot of P vs. B s for the MSPs. The triangles and diamonds represent the field 
isolated MSPs and the field binary MSP with an orbital period smaller than 10 days. In 
addition, the circles shows the MSPs with a characteristic age larger than 10 Gyr. The 
filled square represents the newly born MSP J1023+0038 with B s = 3x 10 8 Gauss, which is 
indicated as an upper limit. The solid, dashed and dotted-dashed line present the relation 
P e = 1.7kB 7 ^ with k — 1,k — 0.4 and k = 0.2. The dotted line represents the line of / = 1 
with E 01 = 1. The results are for i? 6 = M 1A = s± = 1. 
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